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Abstract 

By constructing the commutative operators chain, we derive conditions for 
solving the eigenfunctions of Dirac equation and Schrodinger type equation via 
separation of variables. Detailed calculation shows that, only a few cases can be 
completely reduced into ordinary differential equation system. So the effective 
perturbation or approximation methods for the resolution of the spinor equation 
are necessary, especially for the nonlinear cases. 

PACS numbers: 02.30.1k, 02.30.Jr, 03.65.-w, ll.30.-j 
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1 Introduction 

All fermions have spin-i, which are naturally described by spinors. Many physicists 
such as H. Weyl, W. Heisenberg, once proposed using nonlinear spinor equations to 
establish unified field theory for elementary particlesfl, 2]. Some rigorous solutions 
for the simplest dark nonlinear spinor models were obtained and analyzed [3] -[10]. The 
nonlinear spinor coupling with self-electromagnetic field was calculated in [11]- [15], 
regrettably, only roughly approximate results were obtained due to the complexity of 
the interaction. In contrast, for the linear Schrodinger equation, Pauli equation and 
Dirac equation coupling with external potential V(r) and A(r, 8), there are a lot of 
refined exact eigen solutions obtained via separation of variables[16]-[27]. 
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In the resolution of eigenfunctions to the linear Dirac equation, we find that the 
commutative relations[16] 

[H,J g ]=0, [H,K}=0, [K,J z }=0 (1.1) 

play an important role, where the operators (J Z ,K,H) stand for angular momentum, 
spin-orbit coupling and total energy operators respectively. They are good quantum 
numbers guaranteeing the common eigenfunctions exist. Noticing the definition of 
operators chain 

J g = J z (d v ), K = K(d v ,de), H = H(d v ,d e ,d r ), (1.2) 

which enable us to solve the eigen solutions via separation of variables, and then the 
original problem reduces to ordinary differential equations. In [15], we derived and 
simplified the complete nonlinear Dirac equation system, but where it is still a partial 
differential system depending on coordinate (r,9). The primary motivation of this 
paper is looking for the conditions to reduce the nonlinear Dirac equation into the 
ordinary differential one, and clarify to how much degree we can expect the method of 
separating variables to solve the nonlinear problems. 

Denote the Minkowski metric by rj^ = diag[l, —1, —1, —1], Pauli matrices by 

-MC :)<(: :H: :)} « 

Define 4x4 Hermitian matrices as follows 

'-{(l.'X: :)}-(-)•'•(-)• » 

In this paper, we adopt the Hermitian matrices (1.4) instead of Dirac matrices 7 M , 
because this form is more convenient for calculation. 

Our basic problem is to examine a nonlinear spinor field <fi moving in 4-vector 
potential A^. The corresponding Lagrangian is generally given by[12, 15] 

£ = 0+[ a "(/wfy - eAy) - ^c 7 ]0 + F( 7 , $), (1.5) 

where fi > is a constant mass, F is the nonlinear coupling potential, which is the 
polynomial of 7 and (3 defined by 

7 = + 7 0, P = (1.6) 

It is well known that 7 is a true-scalar and (3 a pseudo-scalar. The variation of (1.5) 
with respect to <p + gives the dynamic equation 

hid (p = H(f), H = a-{-hW -e^ + eAo + ific-F^-f-Fpp, (1.7) 
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where F y = 2£,F« = %. 

Let coordinate x 3 = z along the direction of magnetic field B, then we locally 
have[8, 9, 12, 15] 

A — A(r,9)(— sin ip, cos tp, 0), (1.8) 

— * 

which satisfies the Coulomb gauge V ■ A = 0. In the spherical coordinate system 
(r, 0, </?), we have 

1 1 

<7 • V = ov<9 r + -cr e d e H — <7<o9<o, (1.9) 

r rsinti 

where (cr r ,ag,a v ) is given by 



cos# sinfor^M -sin 9 cos9e~' pi \ / -ie - ^ 
sinOe^ - cos / ' \cos 6e<<* sin ) \ie^ 

Let J be the angular momentum operator for the spinor field 



(1.10) 



J = fx p+-hS, p=-MV, S = diag{a,a). (1.11) 

Then the angular momentum J 3 = J z commutates with the nonlinear Hamilton oper- 
ator (1.7), and the eigenf unctions of J z = —hid v + \hS% are equivalent to the following 
form 

K(pi — ti J (1-12) 

with (k = 0, ±1, ±2, • • • ), where Uk,Vk(k = 1,2) are real functions of (r,9) but inde- 
pendent on if and t, the index T stands for transposed matrix. 
Making variable transformation [15] 

u = u 1 (r,9)+u 2 (r,9)i, v = v x (r, 9) + v 2 (r, 0)i, (1.13) 

then we have 

&o — \ u \ 2 + I 1 2 , ol = (uv — uv)i, ^ 

7 = \u\ 2 — \v | 2 , $ — uv + uv, 

with a = a(— sin ip, cosip, 0). The eigenfunctions of (1.7) with even parity take the 
form 



= U ^ 2n6 ^ v = E Ke ( - 2n+1)01 , (1.15) 



AT TV 

// 

n=-N ra= -JV 

and the eigenfunctions with odd parity take 



AT AT 

// 

n=-N n=-N 



= U n e (2n+1)e \ v = Vne- 2ne \ (1.16) 
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where all (U n , V n ) are real, N < oo for some cases, but N = +00 for general nonlinear 
cases. In what follows we only consider (1.15), because for (1.16) we have similar 
results. By (1.14) and (1.15) we have 

N 

«o = (UnU m + V n V m )cos(2(n-m)6), (1.17) 

n,m=—N 
N 

7 = ^ (U n U m -V n V m )cos(2(n-m)9), (1.18) 

n,m=—N 
N 

a = 2U n V m sm[(2(n + m)-l)0], (1.19) 

n,m=—N 
N 

$ = 2U n V m cos[(2(n + m)-l)6]. (1.20) 

n,m=—N 

Obviously, only for the case N — 0, the quantities d and 7 are independent of the 
following analysis shows only this case can be solved by separation of variables. 

2 Operators and Commutative Relations 

In spherical coordinate system, by straightforward calculation, we have the following 
explicit operator relations[15, 16] 

1 1 

d x = sin 6> cos y9<9 r H — cos 9 cos (pd e -sin</?cL, (2.1) 

r r sin 9 

d v = sin 9 sin ipd r H — cos 9 sin H -cos</?cL, (2.2) 

r rsmfl 

d 2 = cos9d r ^-^d v , (2.3) 

rsm9 

For orbit angular momentum operator L = r x p, we have 

L x = /w (sin + cot 9 cos y%^) , (2.4) 

L y = hi (— cos (pd e + cot 6* sin ^c^) , (2.5) 

L 2 = -hid v , (2.6) 

L 2 = -n 2 (9 2 + cot^ + -^-^). (2.7) 

sin t/ z ^ 
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For the nonlinear system (1.7), detailed calculation shows that 



h 



h 



[H, J x ] = sin tpde (eV - F 7 7 - Fp(3) + e (cot cos 2 ipA + sin 2 ipd e A) a 1 + 

e sin ip cos y?(cot OA — dgA)a 2 — eA cos >pa 3 , 
[H, J y ] = — cos (pdg (eV — F 7 7 — Fp0) + e sin p> cos y?(cot OA — dgAja 1 + 



:2. 8 ) 



h 



[H,J Z 



e(cot sin 2 (pA + cos 2 <pdeA)a 2 — eA sin y?a 3 , 



0. 



(2.9) 
(2.10) 



(2.10) shows J z is still a good quantum number for the nonlinear Dirac equation. By 
(2.8)-(2.10) and relation [H, J 2 } — [H, J] • J + J • [H, J], we have 

[H, J 2 } = h 2 (L 7 + Lp- eL v - eL^simpa 1 - cos ipa 2 ) + D v + eD A ) , (2.11) 

in which 



L 7 = -<9 2 F 7 - cot#<%F 7 - 2d e F 1 d e , 

Lj3 = -d 2 Fp- cot ed e Fp -2deFf.de, 

L v = -d 2 e V -cotOdeV -2d e Vd e , 

L A = -dl A - cot Ode A - 2d e Ad e + cot 2 OA, 

D v = — i(ed e V — <9 e F 7 7) (simp Si — cos (pS 2 ), 



(2.12) 
(2.13) 
(2.14) 
(2.15) 
(2.16) 



D A = (d e A + cotOA)-f(3 + 2AcotO(cotO(cospa 1 + simpa 2 )-a 3 )d<p. (2.17) 

Define the spin-orbit coupling operator by 

K = -f(L-S + h) = h-f- hi(K e d e + K v d v ), (2.18) 

where Kg and fC^ defined by 

-ie-** \ ( ie - "* 
/ ' I -ie^ 



diag 



then we have 



diag 



1 - cot Be"** 
-cotOe^ -1 



-1 cotOe-^ 
cot Oe^ 1 



(2.19) 
(2.20) 



K 2 = J 2 + -K 2 . 
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Therefore, K is the first order form of J 2 . For (1.7), we can check 

[H,K] =K V + K A + K f3 , 



(2.21) 



(2.22) 



in which 



K v = —Hi(edgV — dgF J , y)^(sm ipSi — cos ipS 2 ), (2.23) 

K A = -he (A-f(sm (pa 1 - cos (pa 2 ) - (d e A + cot OA + 2Ad B )(5) , (2.24) 
Kp = 2h-f(3F /3 + h(d e Fp + 2Fpd e ) (sin (pa 1 - cos (pa 2 ) + 

2HF P (cot #(cos (pa 1 + sin (pa 2 ) - a 3 ) d v . (2.25) 

In some sense, (2.22)-(2.25) reflect the influence of parameters on the integrability. 



3 Conditions for Separation of Variables 



For the nonlinear Hamiltonian (1.7), the commutative relation (2.10) shows J z is still 
a good quantum number, but (2.22) shows K is not. In what follows, we look for 
a new operator T and the conditions, such that the systems (1.7) can be solved by 
separation of variables. This is equivalent to the existence of operator T = T{de, d^) 
or T" = T'(d r , dy,), which satisfies commutative relations 



[J 2 ,f]=0, [H,f]=0. 
We only consider the following first order operator 

f = T o -m(T d e + T v d v ), 



(3.1) 



(3.2) 



where {Tq^Tq^T^) are all Hermitian matrices, and their components are smooth func- 
tions of (r,9,tp). The higher order operators can be similarly constructed, but they 
seem not to be functionally independent of (3.2), and the highest order derivatives 
should matches that in H. For Dirac equation, it is enough to look for the linear oper- 
ators similar to (3.2). But for the Schrodinger and Pauli equations, we should consider 
the second order operator similar to L 2 as shown below. Different from p and J, the 
operators constructed from the procedure may have not manifest physical meanings. 
According condition [J z ,f] = 0, we can solve 



T n = 



P 21 cot Be' 1 
P31 

y P 4 i cot 6e l 



Pi 2 cot 9e~ 

P22 
P32 cot 9e~ 

Pt2 



Pl3 

P 23 cot Qe l 

-P33 
P 43 cot Oe' 1 



P14 cot 9e 

P24 
P 34 cot 6e-^ 
P44 



(3.3) 
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/ 



T 



\ 



Ta 



\ 

N u 
N 21 cot Be 1 

n 41 cot 



M n 


M 12 e - ^ 


M13 


M 14 e" J 


M 21 e^ 


M 22 


M 23 e^ 


M 24 


M 3 i 




M 33 




M 41 e i<?i 


M 42 


M 43 e^ 


M 44 



j0 \ 



(3.4) 



N 12 cot 0e -i * 



AT. 



22 



iV 32 cot Be'** 



42 



iVl3 

N 23 cot 6e l 

N33 
N i3 cot 9e i 



Nucotee-** \ 



N- 



24 



iV 34 cot 0e -i * 



N, 



44 



(3.5) 



/ 



where (Pm, Mm, Nm) are Hermitian matrices, and their components are functions of 
(r,6), the factor 'cot#' is introduced for convenience of the following resolution. 
The relation [H, T] can be expressed as 



[H,f\ = H ee d 2 e +H re d rd + (H e + H e ^)d e + (H r + H r ^)d r 
+ (Hq + H^dp + i? w <9 2 ) , 



(3.6) 



where the coefficient matrices can be obtained by straightforward calculation, and they 
can be easily generated by computer. (H 9e , H rd , H 9lfi , H rtp , H w ) are only functions of 
{Pki, Mkh N M ), but the others depend on their first order derivatives. The operator d v 
should be replaced by matrix D v due to solution (1.12), 



D v = i diag(/t, k + 1, k, k + 1), 



d% <-> L>J = -diag(« 2 ,(« + l) 2 ,« 2 ,(« + l) 2 ) 
Then [H, T] = is equivalent to the following equations 



(3.7) 
(3.8) 



Hee — H r 


= 0, 


He + Hq^D^ 


= 0, 




= 0, 




= 0. 



(3.9) 
(3.10) 
(3.11) 
(3.12) 

To solve equations (3.9)-(3.12), we should solve the simple and algebraic relations 
at first, then the complex and differential equations, otherwise the calculation will 
become terribly complicated. The differential equations including in (3.9)-(3.12) are 
linear, and one component usually satisfies two independent equations, which lead to 
constant solutions. So it is easy to resolution. 
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By (3.9), we can solve 



Tg = 



( Wi 

iW 2 e i(/) 
W 3 



W 3 
-W A e* 

W x 
-iW 2 e l 



W A e-' 
W 3 

iW 2 e~ 
W 1 



\ 



(3.13) 



/ 



-iW 2 e~ i4 

W 1 
-Wte-** 
W 3 

in which W k = W k (r,6,K), (k = 1,2,3,4). Again by (3.10), we find that W 1 = W 3 = 
W4 = 0, and W 2 7^ is a constant. Obviously T can be only determined by a constant 
factor and constant translation, so we choose W 2 = 1 to fix the solution. Then in 
equivalent sense we have 

T e = K e . (3.14) 

For the other equations in (3.10)-(3.12), detailed calculation shows that the solu- 
tions for P k[ and N ki are highly underdetermined. Expressing P k[ by N kh we can get 
the following results via long but simple calculation 



V = V(r), A = r- 1 U(9), d B F 1 



0. 



0, 



T - HiT^Dp = (K - hiK^D^) + eC/diag(<70, -a e ) + A , 



(3.15) 
(3.16) 



where A is a constant. If we set f = K when U — 0, we get A = 0. Substituting 

D u 



up > 8^ into (3.16), in equivalent sense, we finally get 
T = K + eC/diag((r e , -a e ), T e = K e , 



T = K 



(3.17) 



In the spherical coordinate system, the above procedure is invertible, so equations 
(3.15) actually form the sufficient and necessary conditions of the separation of variables 
to spinor equation with axisymmetry. Noticing the relation (3.17), conditions (3.15) 
can also be verified by the commutative relations (2.22)-(2.25). 

{V = A = Fp = dgF-y = 0} is just the simplest nonlinear case discussed in [3] -[9]. 

{V = V(r), A = Fp = 0, T 7 = G(r)} are the cases discussed in [16] and [17]-[26]. 



{V 



Ze 



A 



(a — bcos9), Fp = F 1 = 0} is the case calculated in [27]. 



For the general case of (3.15), although the separation of variables is valid for the 
Dirac equation (1.7), but the solutions to the reduced ordinary differential equations 
usually can not be expressed by elementary functions in finite form. So some auxiliary 
numerical computation is still necessary. 

For the Schrodinger equation 



hid (J) = H<j), H = — 



2m 



d r + ^ + ^~ 2 L 2 

r nrr 1 



+ V(r,9), 



(3.18) 



the calculation is simple. Obviously, for L z = —hid^, we have 

[H, d v ]=0, d v ^Ki. (3.19) 
We look for the following second order operator 

Y = Y + Y e d e + L 2 . (3.20) 



By [Y, d v \ = 0, we get 



By [H, Y\ = 0, we find 



d v Y = d v Y e = 0. (3.21) 



V = W(r) + ±U(6), (3.22) 
Y = \ -2mU(6), Y e = 0, (3.23) 

where U{9) and W(r) are arbitrary functions. (3.22) is the condition to solve the 
eigenfunctions of (3.18) via separating variables. Different from the Dirac equation, 
(3.22) shows V can vary with 6. 



4 Some Remarks 

From the above calculation, we get the following insights and experiences: 

(Rl) The above procedure has general significance for separation of variables. The 
method is also suitable for higher dimensional cases. For Hamiltonian 

H n = H n (d 1 ,d 2 ,--- ,d n ), (4.1) 

the separation of variables is equivalent to the existence of the following Hermitian 
operators chain 

H 1 = H 1 (d 1 ), H 2 = H 2 (d 1 ,d 2 ), H n = H n (d 1 ,d 2 ,--- ,d n ), (4.2) 

which form a Abelian Lie algebra 

[Hj,H k ]=0, (j,k = l,--- ,n). (4.3) 

The operators chain {Hk} reflects the intrinsic symmetry of the system, which can be 
obtained via simple calculation. In general, there may not be a complete chain similarly 
to the case of (1.7), but the incomplete chain is still helpful to simplify the dynamic 
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equation as done in [15]. In addition, the existence of operators chain depends on the 
coordinate system. 

(R2) The commutative relations (4.3) are conditional relations, that is, it is enough 
for them to hold only for the eigenfunctions, rather than for the whole function space 
of the solutions. The solution to the dark spinors is an expample[3]-[9]. So we can use 
the information of solutions for solving the operators chain. 

(R3) For the spinor with self-electromagnetic interaction, we have dgV ^ and 
rA 7^ eU(9). Then according to (3.15), we can not construct an operator T commuting 
with H. This means (1.7) has not complete operators chain in the spherical coordinate 
system, so the nonlinear Dirac equation can not be generally reduced to ordinary 
differential equations. Consequently, the equation derived in [16] can not be further 
simplified. Some physical parameters such as g/2 = 1.001159652 •• • may only be 
computed by perturbation. 

(R4) (2.22) provides a measurement for integrability in some sense. By (2.25), we 
find Kp influences many terms, so $ is a very bad parameter for the existence of eigen 
solutions. The dynamical reason for the Pauli's exclusion principle is one of the most 
fundamental problems unsettled. After all, the reason should be something rooting 
in the complete Dirac equation, so the suspicion that (3 is related to this phenomena 
might be reasonable. 

(R5) The above method is related to the concept 'superintegrability', which was 
introduced at the classical level by Wojciechowski[28] and at the quantum level by 
Kuznetsov[29]. Superintegrability also deals with the symmetry and integrability of a 
Hamiltonian system, and many elaborate models were solved. Besides (4.2), the super- 
integrability demands extra commutative operators. There are extensive literatures, 
e.g. [28] -[40] and references therein. However, most of the papers studied classical 
Hamiltonian system and Schrodinger type equations, and few works dealt with spinor 
equations. The above treatment seems to be more straightforward and simpler, the 
information of the solutions and symmetries of the system is automatically used, and 
the resolution seems unnecessarily to introduce extra commutative relations to (4.2). 
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